Tunneling conductance due to discrete spectrum of Andreev states 
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We study tunneling spectroscopy of discrete subgap Andreev states in a superconducting system. 
If the tunneling coupling is weak, individual levels are resolved and the conductance G(V) at small 
temperatures is composed of a set of resonant Lorentz peaks which cannot be described within 
perturbation theory over tunnelling strength. We establish a general formula for their widths and 
heights and show that the width of any peak scales as normal-state tunnel conductance, while its 
height is < 2e 2 /h and depends only on contact geometry and the spatial profile of the resonant 
Andreev level. We also establish an exact formula for the single-channel conductance that takes 
the whole Andreev spectrum into account. We use it to study the interference of Andreev reflec- 
tion processes through different levels. The effect is most pronounced at low voltages, where an 
Andreev level at Ej and its conjugate at — Ej interfere destructively. This interference leads to 
the quantization of the zero-bias conductance: G(0) equals 2e 2 /h (or 0) if there is (there is not) a 
Majorana fermion in the spectrum, in agreement with previous results from 5*-matrix theory. We 
also study G(eV > 0) for a system with a pair of almost decoupled Majorana fermions with splitting 
Eo and show that at lowest Eo there is a zero-bias Lorentz peak of width W as expected for a single 
Majorana fermion (a topological NS-junction) with a narrow dip of width E 2 /W at zero bias, which 
ensures G(0) = (the NS-junction remains trivial on a very small energy scale). As the coupling W 
gets stronger, the dip becomes narrower, which can be understood as enhanced decoupling of the 
remote Majorana fermion. Then the zero-bias dip requires extremely low temperatures T < E 2 /W 
to be observed. 
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der the conditions of these experiments, the Majorana 
state is localized and thus is expected to produce a con- 
ductance peak proportional to the peak in the density of 
states, in agreement with Eq.([T]). For the vortex core, 
direct STM measurements have been propose d 10 ' 16 with 
a similar prediction of a zero-bias conductance peak. 

The idea of probing the Majorana bound state with 
tunneling transport measurements motivated us to re- 
visit the conventional relation ([TJ for the case of a Ma- 
jorana level and, more generally, for an arbitrary dis- 
crete spectrum of Andreev bound states. In this work 
we study the conductance through a tunneling contact 
between a normal metal and a superconductor hosting a 
This relation allows to probe local density of states with subgap spectrum of Andreev bound states Ej . While the 



I. INTRODUCTION 

Scanning Tunneling Microscopy (STM) is a powerful 
method of studying electronic surface properties of con- 
densed matter systems. Electrons from the metallic tip 
of the microscope tunnel into the system, creating a sta- 
tionary current at fixed voltage. At low temperatures the 
conductance G of the tunnelling contact is proportional 
to the surface density of states p beneath the tip, and the 
coupling strength between the tip and the system 



G{V) cx \t 2 \p(r,eV). 
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energy and spatial resolution. In recent works*— trans- 
port measurements have been suggested to probe the 
Majorana fermion in superconducting structures with 
non-trivial topological properties, such as nanowires 
combining strong spin-orbit, induced superconductiv- 
ity and magnetic field"^, and vortex cores on 
topological insulator surfaces with induced s-wave 
superconductivity- 7 - ' 10 ' 15 ' 16 . There are two major trans- 
port phenomena suggested, the first one is the fractional 
Josephson effect* 1 -' 7 - ' 14 ' 17 : the current-phase relation in 
an SNS-j unction hosting a pair of Majorana fermions is 
4-7r-periodic with phase if the fermionic parity in the junc- 
tion is conserved 1 ^. This effect also implies the doubling 
of Shapiro steps in an ac measurement^ 1 -, which has re- 
cently been reported 1 ^. A second, and more direct ef- 
fect of the Majorana state on transport is the peak in 
zero-bias conductance. Transport experiments have been 
done on nanowire-based setups^ - - — , showing a zero-bias 
peak, that is robust to the change of magnetic field (in 
both magnitude and direction) and gate voltages. Un- 



G{V) characteristic is expected to reflect the structure 
of the spectrum and exhibit peaks at voltages close to 
Ej/e, it is clear, that the classical expression (JlJ cannot 
accurately describe the current in our system of inter- 
est. Since we consider a discrete spectrum, the density 
of states p(E) is a set of (5-functions. The conductance 
however, cannot be a set of 5-functions as it is limited 
by the number of channels in the contact. Moreover, ex- 
pression ([l} is a Fermi Golden rule result, that describes 
electrons slowly leaking the probed system. In our su- 
perconducting system there are no extended states at 
subgap energies and single electrons cannot leak into it. 
The only possibility for a stationary current is due to An- 
dreev reflection: the process of electrons being reflected 
from the NS-junction as holes. The discrete Andreev lev- 
els in the superconducting system act as quasi-stationary 
levels for this process: an incident electron may occupy 
an energetically favorable level and is then emitted back 
into the metal tip either as an electron or as a hole. At 
voltages close to one of the levels Ej this process should 



2 



exhibit resonant behavior, and cannot be described per- 
turbatively in terms of the coupling strength. Our aim 
is to describe these resonances for a generic tunneling 
contact. 

Another interesting problem of the discrete spectrum 
is interference of Andreev reflection processes involving 
different levels. There is no interference involved in the 
simple formula (TTJ) , where the conductance is directly pro- 
portional to the number of available energy levels on the 
probed surface. The process of Andreev reflection is dif- 
ferent from conventional tunneling as it involves tunnel- 
ing into the probed system and then back into the metal. 
Reflection of the electron into a certain hole channel can 
happen in multiple ways, that differ in the Andreev level 
Ej they tunnel into. The amplitudes of these different 
ways have to be summed together, which allows for in- 
terference between different levels. The interference ef- 
fects should be important for a few-channel contact, es- 
pecially a single-channel one, where all amplitudes are 
added together. Most interesting is the case of low volt- 
ages, where the BdG-symmetry of the spectrum becomes 
important and defines the interference between conjugate 
levels ±Ej. 

In section|TT]we formally introduce the system we study 
and the formalism we will use. This includes the Lan- 
dauer formula for superconductors, a Kubo-type formula 
for Andreev reflection, and the T-matrix expression of ex- 
act Green's functions in terms of the tunneling Hamilto- 
nian and unperturbed Green's functions of the connected 
systems. In section lTni we study the resonance peaks that 
occur at voltages close to a single Andreev level Ej . We 
derive general formulae describing the peaks and study 
the particular cases of a peak on a Majorana level, for- 
mulae for a point-contact and the conductance at finite 
temperatures. In section IIVI we establish an exact for- 
mula for conductance in the single-channel contact case 
that takes the whole spectrum Ej into account. We then 
study interference of different levels and observe quan- 
tization of zero-bias conductance. Next, we study con- 
ductance at finite voltages. In particular, we consider 
a system with a pair of weakly coupled Majorana lev- 
els. We find that a zero-bias peak in this system coexists 
with an extremely narrow dip at E = that shows that 
the system is still topologically trivial. We discuss the 
experimental requirements to observe the dip structure 
of this peak and show that, surprisingly, they become 
more severe when the tunnelling coupling to the probe 
increases. In section [V] we discuss some experimentally 
realized systems where level spacing of Andreev states is 
large enough for the individual levels to be observed. 



II. SETUP 



electrons to tunnel from the tip onto the system. 

H = H M + H T + H s (2) 

We are interested in voltages and temperatures much 
lower than the superconducting gap and thus only take 
the Andreev levels into account: 



(3) 



The contact is described by the tunnelling Hamiltonian 
Ht = 7" \0 m )t ms (a s | + h.c. (4) 



where \0 m ), \a s ) describe some electronic wave functions 
in the metallic tip and superconductor, correspondingly. 
In the Bogolyubov-de Gennes formalism we will use, Uj 1 
has to be extended into Nambu space 
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where is the time-reversal operator, and |^) 
is the particle-hole-conjugated partner of \0 m ), i-c. 
\9* n ) — C\9 m ) with C denoting particle-hole conjugation. 
Throughout the paper the check sign indicates a struc- 
ture in Nambu space. 

At subgap energies E < A electrons incident from the 
metallic tip are reflected back either normally or as holes. 
The latter process is Andreev reflection and produces a 
current through the NS-junction, which can be described 
in the Landauer formalism. The formula describing the 
stationary current at some voltage bias reads (we derive it 
from the general Lesovik-Levitov formula^ for the BdG- 
Hamiltonian) 

oo 

I=j J T A (E)[f (E - eV) - f Q (E + eV)]dE (6) 
o 

Note that the factor of 2e is due to the quantum of trans- 
ferred charge in the Andreev process. fo(E) denotes the 
equilibrium Fermi distribution function, and T A (E) is the 
sum of Andreev reflection eigenvalues: 



T A (E) = tr (S eh (E)^S eh (E)) 



(7) 



where S e h(E) is the electron-hole block of the ^-matrix of 
the N-S junction. With the help of the relations Ta(E) = 
T A (-E) and f (E) + fo{-E) = 1 formula © can be 
rewritten in a form similar to the one for the standard 
N-N contact current: 



oo 

/ = X/ T A(E)[fo(E - eV) - f Q (E)]dE 



(8) 



The system we are studying consists of a metallic tip 
described by the Hamiltonian Hm, a superconducting 
system Hs and a tunneling Hamiltonian Ht that allows 



dl_ _ 2f_ 
~ dV ~ h 
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Below we mainly study the zero-temperature conduc- 
tance G(eV) — T A (E — eV) ■ 2e 2 /h which is just the sum 
of Andreev transmission eigenvalues times the relevant 
conductance quantum 2e 2 /h. Conductance at T > can 
then be easily obtained from the zero-temperature result 
via the convolution (jSJ). The trace (0 can be written as 



T A = tr 



v 




Gi 





v* 



G 



(10) 



which can be understood as an Andreev-reflection ana- 
log of the Kubo formula for the linear response, v is 
the velocity operator for electrons and v* — 9w9 _1 is 
its time-reversal. The explicit matrices inside the trace 

_____ v A R 

([TO]) act in Nambu space, and G E are exact advanced 
and retarded Green's functions at energy E. The exact 
Green's functions in the metal are exressed through un- 
perturbed Green's functions G = [E + iO — H M ]~ X and 
Gs — [E + iO — Hs] 1 (here and in what follows Green's 
functions without index are assumed retarded): 



T 



G — Gq ~\~ TGq 
Gq Hx Gs Ht 



(ii) 



fl2) 



1 — GqHtGsEIt 
The Green's functions here act in the Nambu space, so 

'E + iQ-Hu V 1 (g E 
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E + iO + H 



M 



9E 

-9-E 
(13) 



We use g to denote normal metal Green's functions with- 
out Nambu structure and g* — QgQ^ 1 for its time- 
reversal image. The superconducting Green's function 
is 



liXil 



j>0 



E + iO-Ej E + i0 + Ej 



(14) 



The sum over j > reflects the symmetry of the levels 
- each positive level \j) has a particle-hole conjugated 
partner \j*) with opposite energy. Gs has poles at E = 
±Ej and the T-matrix cannot be treated perturbatively 
in their vicinity. All terms of the geometrical sequence 
have to be included in the expression (|12[). This can be 
done in a compact way in two cases we consider below. 



III. SINGLE-LEVEL RESONANCE FOR 
ARBITRARY H T 

The T-matrix and consequently Ta can be found 
explicitly in a closed form in the resonance situation 
\E — Ej \ -c \E — Ej^j \ when the energy of the inci- 
dent electron is closer to the Andreev level Ej Q than to 
any other level Ej. We only keep the (E — Ej ) _1 term 



in Gs in this case, and write 



HtG sHt 



G \t)(t\ 



E - Ej, 



iO 



with \t) = H T \jo) 



(15) 



Here |r) is a function of the electron coordinates in the 
metallic tip, which is determined by the contact Hamil- 
tonian Ht and the resonant Andreev state |jo). It has 
both electron and hole components which we denote as 
r e ) and \rh) for later use. The simple projector form of 
the operator (p~5)l leads to 



T = 



Go\t)(t\ 



E-E jo +iO-(r\G \T) 



(16) 



The real part of (t\Go\t) is zero, and we denote its imag- 
inary part as W. Substituting the T-matrix (|16|) into 
formula (fTTj) and the result into the Kubo formula (|10[) . 
we find 



2e 2 
G=^-T A 
h 



\g E vg^\r e )(T h \g%v*g^ E \T h ) 



(E-E jo y + W 2 



W — n(n e + nh) with 
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(17) 



(18) 



(here |r^) = Q\rh) is just the time-reversed wave func- 
tion) In the case of a clean metallic tip (and assuming 
that Fermi energy Ep is the largest energy scale), the 
gvg-terms can be simplified in terms of the n ei h densi- 
ties: 



G(eV) 



2e 2 



^ rri* 

1 H w 2 — 

_ An e n h 
A — 7 ; ^ 



(19) 



(20) 



Formula (|19[) constitutes our main result for the single- 
level resonance, and we will now discuss it in more de- 
tail. The ^-dependence of (|T9]) has a Lorentz shape, 
which is generic: at resonance, the 5-matrix is governed 
by the pole at E — Ej + iW. The width W describes 
the broadening of the Andreev level at the energy Ej 
due to coupling to the tunneling probe. T* A < 1 is the 
probability of Andreev reflection for an incident electron 
with E = Ej . It is defined by the relation between 
the electron and hole components of \jo) and the geom- 
etry of the contact. If | jo) was a purely electronic state 
T A ~ n h = an d there would be no Andreev reflection 
from | jo ) and hence no conductance. For deep subgap 
states (Eq A) electron and hole components are very 
similar (e.g. connected by phase rotation for a uniform 
A(r)) so we expect 1 — T\ <C 1. However, such a feature 
may be absent when measuring tunneling conductance at 
points where the gap A is locally suppressed, A(r) < Eq, 
e.g. near the origin of a vortex core, or in the middle of 
an SNS-junction. In such regions electron and hole com- 
ponents may differ considerably, leading to a small T\. 
The presence of the T\ factor in (fl"9)) is an important 
difference from the classical relation ([T]) : Andreev reflec- 
tion amplitude depends on the e — h structure of the wave 
function, and not just on the overall probability density. 
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A. Majorana peak 

An important special case of the single-level resonance 
described by (fH?|) is the resonance on a Majorana bound 
state. In this case \jo) is self-conjugate, so that |r e ) = 



Hence n e 



rih and TJ 
2e 2 



1. 



G 



1 



(21) 



This agrees with earlier derivations of the Majorana- 
induced zero-bias tunneling conductance peak^^. The 
width W of the peak still depends on contact and wave 
function geometry, while the height is fixed to 2e 2 /h by 
the symmetry of the Majorana state. Other energy lev- 
el 8 \j jo) generally lead to a negative correction to the 
conductance of the order of 0(W 2 /w 2 ) where luq is the 
typical level spacing of the Andreev spectrum in Eq. (J3j . 



B. Point-contact 

The case most relevant to a scanning tunneling mi- 
croscope setup is the point contact Ht = \6)to((j\ + h.c. 
where wave functions \9) and |cr) are localized at some 
points rjf and r$ in the metallic tip and superconductor 
correspondingly. For simplicity we assume here that Ht 
commutes with spin, but our results are easily general- 
ized to a spin-sensitive setup. 



T* A 



w 



= 7r|t 



(22) 
(23) 



Here vm is the density of states in the metallic tip per 
spin projection and \^ e {h) 2 \ = \ip 2 {h)t \ + \^\ h)i \ is the 
probability density of the electron (hole) component of 
the resonant Andreev level. We can rewrite ([21]) for the 
point-contact as 

.,2 



G(eV) = AK 2 \ti\-v M i>T* A , 



(24) 



which has the form of a standard expression for normal 
conductance times the resonant reflection probability T\. 
v in Eq. (l24l) is a broadened density of states of the An- 
dreev level | jo) : 



l + (eF 



The point-contact expression 
an arbitrary contact: 



G = -|^tr (^Ht [Gq — Gq] Ht 



-E j0 )*/W* 
can be generalized to 



Gq 



T* A (26) 



where Gs is the broadened Green's function of the probed 
system: 

1/(ttW) 



G(eV) = 



1 + (e - eV) 2 /W 2 



G(e)de (27) 



If instead of Gs in (|26|i we substitute the Green's func- 
tion corresponding to the normal state of our probed 
system, and put T* A = 1, the expression ([24| will yield 
the normal-state conductance. In the superconducting 
regime in presence of the Andreev spectrum, the Green's 
function pole has to be broadened according to Eq. (f27|) , 
so that (5-peaks in G R — G A are smeared into Lorentz 
peaks of the width W. We remind that the parameter 
W is different for different peaks according to Eq. (|18[) . 
Thus, if W is calculated for |jo), then expression (|27|) 
with that W only works for the peak at eV = Ej . 



C. Temperature dependence 

Given the zero-temperature G(T = 0,eV) from (Q3 
the finite-temperature conductance is readily obtained 
via formula ([9]) by convoluting it with the thermal distri- 
bution function. The shape of G(T) starts to differ from 
the G(T = 0) case as soon as T becomes comparable to 
W. At this point the peaks become broader, while re- 
taining spectral weight. At luo ^> T 3> W (we denote the 
level spacing by luq) the temperature-broadened resonant 
peak looks like 



2e 2 

G(eV) = —TX 



nW 



Tcosh 2 



Ej -eV 
2T 



(28) 

Here W enters as the spectral weight of the original 
Lorentz peak ([TO]) , while the peak shape is defined by 
the temperature T only. If instead of expression (fl9|) 
we use expression (f24| or (|26|) in the convolution ((9]), we 
find that the broadening of the <5-peaks by W is much 
less than thermal broadening and thus drops out. For 
the point-contact 



G(eV) = TXn 2 



\t>m/2 



h Tcosh 2 (i 



2T 



(29) 



IV. SINGLE-CHANNEL CONTACT, 
MULTI-LEVEL SYSTEM 

Another case which can be solved exactly is presented 
by a single-channel contact Ht coupled to an arbitrary 
system of Andreev levels H$ ■ In this section we account 
for all Andreev levels simultaneously in an exact way, 
which allows us to study important interference effects 
due to Andreev reflection via different levels of the dis- 
crete spectrum. 

An arbitrary single-channel tunneling Hamiltonian can 
be written as \it = \9)t(a\+h.c. with a real t and normal- 
ized electronic wave functions \9) and \<j) that are defined 
in the metal and superconductor respectively. To include 
Ht into formula (|12[) . we extend it into Nambu space: 



H T = \9}t(a\ - \9*}t* (a*\ + h.c. 



(30) 
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Thus, in the BdG formalism Ht appears as a two-channel 
contact. The calculation of the T- matrix is similar to 
that of the expression (fT6| . Using (|30|) and (Q~4|) we write 



G H t G s Ht = G (\0) \6*)) N^t 



with a 2 x 2 matrix E 



E = 



Ea 



y -t 2 N vlifMJl 



E + E, 



E 



h 



IW>I 2 , IHi)l 5 



-t 2 N, 



£ + E, 



i>o 



E + Ei 



(31) 

(32) 
(33) 

(34) 

(35) 



Here Ao = lm(8\g§\8) is a normalization factor (for a 
point-contact Nq — ~nv) . E^ has a simple physical mean- 
ing: 2E^ is the amplitude of Andreev reflection (up to 
some phase factor) in the lowest order in Ht- Simi- 
larly, E e) /j describe the corrections to electron-electron 
and hole-hole reflection due to the presence of Ht- Us- 
ing Eq. pip we now construct the T-matrix: 



f = G (\9) \9*)) 



1-it 



(36) 



where the 2x2 matrix in Eg. (|3"61 reads as 



1^1 



E R E> 



1-iS 



(l-iE e )(l-tE h ) + |E= 



(37) 



Finally we substitute G = G + TGo (and G A = G Rjf ) 
with T-matrix (|36|) into the Kubo formula (|T0| . Only 
off-diagonal (in Nambu space) terms of G R ' A and conse- 
quently only off-diagonal terms of (f3"7| enter the trace, 
yielding 



G( e y) 



2e 2 



4|Si| 



(1 - E e E, 



(38) 



Voltage eU enters the formula through E where E is set 
equal to eV. At resonance, |eV — Ej \ <C IU, formula 
(|38|) reproduces our single-level result (|19|) with n e = 
|(cr|j)| 2 t 2 A" /7r and n h = \(a\j*)\ 2 t 2 N /ir. Away from 
any resonance, E <C 1, so that G(eV) ~ 2e 2 /h x 4|E,4| 2 , 
which is the lowest order perturbative result in powers of 
Ht- Ea is a sum of Andreev amplitudes over different 
levels j, and at some voltages several levels contribute to 
T, A terms of comparable magnitudes, which may inter- 
fere. These interference effects are most drastic at lowest 
voltages, where interference pushes G(0) to be equal ei- 
ther to or to 2e 2 /h, as we discuss below. 



A. Zero-bias conduction quantization 

The fact, that the zero-bias conductance of a single- 
channel NS-junction is quantized follows from general 
symmetry consideratio n 24 ' 25 . The ^-matrix describing 
an NS-junction has BdG-symmetry which relates and 
S-e- In conjunction with unitarity S'S = 1 this fixes 
det = ±1, dividing possible S-matrices into two dis- 
connected classes: trivial (det So = 1) and non-trivial 
(det So = — 1) . The non-trivial case is represented by 
topological superconductor-metal junctions supporting a 
Majorana level on the junction 2 -, (this Majorana level is 
a well-defined quasistationary Andreev level if there is 
a barrier in the metal in the vicinity of the contact to 
superconductor, otherwise the Majorana level is leaked 
into the normal metal). In the single-channel case S is a 
2x2 matrix and G(0) = (1 - det S )e 2 /h. Thus G(0) is 
(2e 2 /h) in the trivial (non-trivial) case. 

Let us now show that expression (I38p is indeed quan- 
tized at E = 0. Suppose first there are no zero levels 
among Ej. Then, E^ = since terms with |j) and |j*) 
in (|35|) cancel each other out at E = 0, and consequently, 
G(0) = 0. Suppose now that there is a single Majorana 
zero level jo- Then \j ) = |jg) with E = 0, leading to 
E e ~ E ft ~ |E A | - 1/E and E e E ft - |S^| = 0(1) in the 
limit E — > 0. Consequently, G(0) = 2e 2 /h in this case. 

The two above cases (no Majorana states, or a sin- 
gle one) describe all physical situations, as any higher 
number of zero modes would split. Let us show how it 
happens in somethat more details. One possible case cor- 
responds to a non-Majorana zero level with Eq = and 
| jo) 7^ |jo), which can also be understood as a pair of un- 
coupled Majorana fermions |jo) = e ltp {\li) + £|7a))- De- 
pending on the choice of ip, the state |jo) can be decom- 
posed into different pairs of Majorana fermions 71, 72- If 
IHio>| = IMJo*)I we choose e 2 ^ = (a\j ) / (*\j* ) and 
find that the correspoding 72 is not coupled to the con- 
tact, (a\j2) = 0. Thus, 72 does not enter the hamil- 
tonian H = Ht + Hm + Hs an( l we en d up with the 
case of a single Majorana fcrmion 71. In the second 
case, I (c I jo) I 7^ K^b'o)! both and 72 enter H T for any 
choice of (p. The Majorana fermions are effectively cou- 
pled by the contact into a finite-energy fermionic level, 
which drives our NS-junction into the trivial regime. In- 
deed, if |(cr|jo)| ^ l(cbo)l the 1/E A singularities in the 
denominator of expression (f3"5|) do not cancel out so that 
G(0) = 0, as expected for a trivial junction. Finally, 
one may ask what happens if there are multiple zero lev- 
els jo, jo • • • ■ This problem reduces to those considered 
above by changing the basis of the zero levels so that only 
one level couples to the contact, which is always possible 
for a single-channel contact. 

Note that exact quantization only happens in the 
single-channel case. For instance the exact cancellation 
of |j) and |j*} Andreev amplitudes is only possible when 
the incoming electron and outgoing hole are in the same 
channel. This suggests that there may be some traces 
of the destructive interference for few-channel junctions, 
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which become negligible as the number of channels be- 
comes large. 



B. Conductance at finite voltage and interference 

The quantization discussed above only concerns E = 
eV = and cannot provide any information about finite 
energies. Our formula ([55]) allows to find how narrow the 
quantization region is, i.e. at which voltages (or temper- 
atures) it should be visible. The study of (|38)l at finite 
E also resolves an apparent paradox similar to that dis- 
cussed irt 2 ^ consider a finite superconducting system de- 
signed to host Majorana fermions. e.g* 2 -^. It is topologi- 
cally trivial, because there is an even number of Majorana 
fermions across the system and there is finite energy split- 
ting between them. So, detS = 1 and G(0) = 0. Let 
us now extend the system to move one Majorana level 
to infinity. It is now uncoupled from the system and we 
are left with a topological junction which should show 
a zero-bias conductance peak. However, from continuity 
Sq should remain in the trivial class and G(0) should stay 
at zero. Below we show that both statements are true: 
there is a peak at low energies, and, at the same time, 
G(0) = 0. 

The low- voltage conductance is governed by the lowest 
energy level \j) and so we omit other levels in what 
follows - they only lead to minor corrections in G. We 
find 



G{V) 
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The shift in the resonance energy in (|40l) is negligible, as 
it is smaller than Eq by a factor of W/ujo (except for the 
exotic case mentioned earlier when there are two uncou- 
pled Majorana fermions in the superconducting system, 
which are then hybridized by the tunneling contact; in 
this case Eq is the effective energy splitting and cannot 
be neglected). The two- level formula (|39[) agrees with 
what has been derived earlier in a different way for the 
special case when Ht couples the metal tip to a single 
Majorana operator—. If Eq W (gray curves on Fig. 
11 12 j) . the function (j39| shows two well-separated Lorentz 
peaks at ±E with a width W each. If we decrease E , 
these two peaks move towards each other, but G(0) stays 
zero. When E becomes comparable to W (blue curves 
on Fig. Ill2p the shape does not look like a pair of Lorentz 
peak anymore. In the region \eV\ < Eq the conductance 




Figure 1: (Color online) Single-channel conductance at T = 
for fixed W and various values of En. At En = 5W (gray 
curve) there are two well-separated Lorentz peaks of width 
W each. At En = W (blue curve) the peaks are close to each 
other and significantly deformed. At Eq = 0.15W (red curve) 
the peaks merge into one Lorentz peak of width 2W with a 
narrow dip of the width Eq /W in the middle. At E = 
(black dashed curve) there is no dip, and G(eV) is an exact 
Lorentz peak of width 2W. The inset shows the same picture 
at lowest voltages. TJ is put to unity throughout our figures. 



is suppressed due to destructive interference, while at 
eV > \E \ the conductance is enhanced: at eV » En 
one gets G(eV) = 2e 2 /h ■ AW 2 /(eV) 2 which is twice the 
sum of two Lorentz tails. Finally, at Eq <C W (red curves 
on Fig. I1I2[) . we see one accurate Lorentz peak of dou- 
ble width 2W^ with a parametrically narrow dip in the 
region E < E 2 /W all the way to 0. While quantiza- 
tion G(0) = is fulfilled, it only affects a very narrow 
region of energies. Note that throughout the change of 
En the current I(V > E ,W) = G(V)dV (propor- 
tional to the area under the curves on Fig. [1]) remains the 
same, and equals 2e/h-irW. Hence, at high temperatures 
T ^> Eq , W the interference effect upon the conductance 
is exponentially small. 

The above consideration solves the apparent paradox: 
the region of voltages and temperatures where the triv- 
ial character of the system is observable is very nar- 
row, eV,T < Eq, Interestingly, increasing the coupling 
strength W to the probe makes the observation of the 
conductance dip even more difficult - at W 3> Eq one 
needs to resolve energies of the order of Eq/W instead 



of just Eq, as can be seen from Eq. (131)1) and Figs. I1I2I3"! 
This can be understood in terms of the specific topolog- 
ical transition of the S'-matrix described by Pikulin and 
Nazarov 25 . Our S'-matrix has two poles (which are also 
poles of function at E± = ±y/E^ - W 2 + iW. For 

a weak contact (gray curve on Fig. [2]) they are connected 
by BdG-symmetry, E+ = —E*_ and represent a quasis- 
tationary Andreev level at E+ and its image at £L . As 
W is increased and the coupling gets stronger the poles 
move towards the imaginary axis. At W — Eq (blue 
curve on Fig. [2]) they meet at the imaginary axis and 
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Figure 2: (Color online) Single-channel conductance at T = 
for fixed Eq and various values of W. At W = 0.2_Eo (gray 
curve) there are two well-separated Lorentz peaks of width 
W each. At W = Eo (blue curve) the peaks are close to each 
other and significantly deformed. At W = IQEq (red curve) 
the peaks merge into one Lorentz peak of width 2W with a 
narrow dip of the width Eq/W in the middle. 
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Figure 3: (Color online) Single-channel conductance at Eo = 
1 and W — 10 for different temperatures. At T = 0.01 
(red curve) there is a sharp dip with a finite but small G(0). 
At T = 0.1 (blue curve) which is of the order of the zero- 
temperature dip width Eq/W , the dip is significantly smeared 
out, while the wings of the Lorentz peak remain unchanged. 
At T — 1 (gray curve) which is of the order of Eo there is 
only a very shallow dip left. Finally, at T — 4 (black dashed 
curve), which is comparable to W = 10, there is no visible 
trace of the dip and the Lorentz peak starts to broaden. 

the transition^ 5 , happens. At higher W one of the poles 
moves up the axis, the other down and each of them is 
now BdG-self-conjugate: E± — —E±. At very high W 
(red curve on Fig. [2]) the lower pole is "buried"— - it sits 
parametrically close to the real axis, _E_ = iE^/ZW , and 
is very weakly coupled to the probe. This means, that 
the better the contact, the more our NS-junction looks 
topological from the electronic transport point of view. 
This is an important observation - while the quantiza- 
tion of G(0) happens for any W, a coupling stronger than 



the splitting Eq diminishes the temperature and voltage 
range in which the conductance dip (and hence the re- 
mote second Majorana fermion) is observable, as shown 
in Fig. m 

V. DISCUSSION 

We now discuss real systems where the effects stud- 
ied can be observed. A natural object to look at is a 
core of an Abrikosov vortex in a quasi-2D superconduc- 
tor. The above results are applicable if the Caroli-De 
Gennes-Matricon (CdGM) level o 30 i 31 accessible by the 
tunnel contact on the top plane of the superconductor 
are not spreaded into the continuous band dispersing 
along the vortex core. The CdGM spectrum of a clean 
vortex core is described by 

P 2 

= + ( 44 ) 

with uiq ~ ztj and m* — m z (pF±£,) 2 ; where \x is the 
angular momentum, £ is the coherence length and m z 
is the effective mass along the z-direction assumed to be 
larger than {x, y)-plane effective mass m in an anisotropic 
superconductor. In the presence of even a low concentra- 
tion of impurities, CdGM states in a neighboring atomic 
planes of layered material will be displaced by the amount 
of energy ~ wq, whereas bandwidth of each level is about 
fiujQm/m z . Therefore, in strongly anisotropic materials 
with m z 3> m individual low-lying CdGM levels are well 
localized 32 ' 33 and tunnelling from the STM tip into the 
surface of such a material may show the features dis- 
cussed above in this paper. However, in order to discern 
these features, one needs to work at very low tempera- 
tures T < ujq. Currently, it should be possible to reach 
this condition while dealing with a material like NbSe2, 
whereas older dat a 34 ' 35 were taken at too high tempera- 
tures. 

Another possibility is to use genuine two-dimensional 
electron systems like graphene or a surface of topolgi- 
cal insulator. In these materials there are no transverse 
quantum numbers and the Andreev spectrum is discrete 
with level spacing u>q (in graphene the spectrum is four- 
fold degenerate due to valley and spin degrees of free- 
dom) . An additional advantage of these two-dimensional 
systems is a controllable Ef. The charge neutrality point 
in graphene coincides with the i-L -point of its Dirac spec- 
trum, which allows to have u>o comparable to Ef. An im- 
portant issue is the implementation of superconductivity 
in these two-dimensional systems. A possible solution is 
to cover them with a conventional superconductor to in- 
duce superconductivity by proximity effec t 15 ' 26 " —. How- 
ever, one needs to make sure that there are no parasite 
Andreev states in the covering superconductor. This can 
be done by making a hole in the superconductor of a ra- 
dius of the order of the coherence lengt h 10 ' 16 ' 29 physically 
removing the normal vortex core with all its low-lying 
Andreev states from the superconductor. In addition, 
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such a hole should facilitate the application of a tunnel- 
ing microscope to the graphene or topological insulator 
surface. 

VI. CONCLUSION 

In conclusion, we have studied resonant Andreev re- 
flection in an arbitrary tunneling N-S contact between a 
metal and a superconductor hosting a discrete Andreev 
spectrum Ej . If the level spacing is larger than the in- 
verse dwell time W, the subgap conductance G(V) con- 
sists of Lorentz peaks centered at eV — Ej. We have 
derived formulae describing these peaks in the general 
case. The heights of the peaks do not scale with the 
coupling strength and are generally very close to 2e 2 /h, 
while widths scale as normal conductance of the con- 
tact. We also established an exact formula describing the 
Andreev conductance for a generic single-channel con- 
tact that takes into account all Andreev levels at once. 
This formula reveals the role of interference between An- 
dreev reflection processes involving different levels Ej. 
At eV — interference between the level Ej and its 
particle-hole image — Ej is destructive, so that G(0) is 
exactly 2e 2 or depending on whether there is a Ma- 
jorana fermion in the spectrum, in agreement with S- 
matrix theory of topological NS-junctions. We have also 



studied G(V) for finite energies for different W and low- 
est energy level Eq and explained what happens when the 
system is gradually driven from the trivial to the topo- 
logical phase by removing a Majorana fermion from the 
system and thus driving Eq to zero. We have found that 
at Eq <C W the conductance shows a zero-bias Lorentz 
peak of the usual width 2W, with a very narrow dip 
on top of it, which secures G(0) = 0. The parametri- 
cally small width Eq/W of the dip can be understood 
as the "burying" of the remote Majorana fermion: the 
strong coupling of the tunneling tip to one of the Majo- 
rana states effectively decreases the coupling of the sys- 
tem to the other Majorana state. To observe the dip 
(i.e. the "buried Majorana fermion") experimentally, 
temperature and energy resolution has to be smaller than 
Eq/W. The effects we discussed should be observable 
in two-dimensional electronic systems including graphene 
or the surface of three-dimensional topological insulators 
with proximity-induced superconductivity. Another pos- 
sibility is to study the level structure of Abrikosov vor- 
tices in a layered strongly anisotropic superconductor like 
NbSe2 at ultra-low temperatures. 

We thank Alexey Ioselevich, Pavel Ostrovsky and 
Dmitry Pikulin for very helpful discussions. This work 
was supported by the RFBR grant # 10-02-00554-a and 
by the RAS program "Mesoscopic and disordered sys- 
tems" . 
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